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PREFACE 


This  report  summarizes  the  research  performed  under  Contract  Nonr~710(22) 
for  the  year  ending  31  Hay  1966,  The  work  is  in  various  stages  of  completion 
as  reflected  by  the  differing  format  of  the  four  reports  covering  the  various 
areas  of  research. 

The  first  two  papers  of  this  report  are  essentially  scientific  reports 
representing  the  culmination  of  one  phase  of  the  research.  The  last  two 
papers  are  more  in  the  nature  of  status  reports  on  current  research.  A  final 
note  on  stellar  scintillation  has  been  appended  simply  to  indicate  continuing 
interest  in  this  field  although  the  demands  made  on  experimenters'  time  on 
other  projects  precluded  further  work  on  this  area  during  the  past  year. 
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An  Experimental  Study  of  the  Equation 
of  Motion  in  the  Free  Atmosphere 
by 

Horner  T.  Mantis 

1.  Introduction 

The  energy  of  atmospheric  motions  is  distributed  over  a  sufficiently  broad 
band  of  frequencies  or  scales  that  a  quantitative  application  of  the  conservation 
laws  of  mechanics  to  a  specific  frequency  domain  may  still  present  difficulty. 

The  conventional  formulation  of  the  equation  of  motion  as  applied  to  the  weather 
map,  for  example,  which  may  be  written 

If  ♦  f  k  x  V  =  -g  5  p  h  +  Fr  (1) 

is  based  on  Reynold's  hypothesis  that  it  is  possible  to  define  the  velocity, 
accelerations  and  forces  on  a  uniform  scale  (Miller,  1951).  The  frictional 
force  Fr  in  this  equation  represents  a  transfer  of  momentum  frcrr  weather  map 
or  synoptic  scale  by  motions  of  lesser  scale. 

Difficulties  with  the  application  of  Eq.  (1)  may  be  appreciated  by  noting 
that  in  the  free  atmosphere  it  has  not  been  possible  to  use  the  equation  even 
as  a  phenomenological  definition  of  the  frictional  force  by  a  computation  of 
a  difference  in  the  forces  and  the  accelerations.  The  problem  is  largely 
that  of  smoothing  the  fields  of  wind  and  pressure  on  the  same  scale  and,  in 
practice,  is  caused  almost  in  equal  measure  by  the  ill  defined  error  of 
measurement  and  incomplete  specification  of  the  wind  and  pressure  by  the 
present  network  of  meteorological  observations. 

This  paper  presents  the  results  of  experimental  study  to  determine  the 
representative  magnitudes  of  the  various  terms  of  Eq.  (1)  as  applied  to  motions 
of  synoptic  scale  in  the  free  atmosphere,  specifically  at  the  500  mb  and  300  mb 
levels.  For  the  reasons  detailed  above,  it  was  necessary  to  attack  the  problem 
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obliquely  by  a  statistical  analysis  of  data  derived  by  objective  procedures. 

As  in  most  statistical  studies,  the  representative  magnitude  of  a  quantity 
is  assumed  to  be  specified  by  its  mean  and  variance.  In  this  sense,  it  has 
been  possible  to  estimate  the  magnitude  of  the  following  quantities  of 
meteorological  interest: 

a)  The  atmospheric  accelerations  at  synoptic  scale. 

b)  The  frictional  force  at  synoptic  scale  (actually,  evidence 
is  presented  simply  to  the  effect  that  the  frictional  force 
must  have  the  sane  magnitude  as  the  synoptic  acceleration). 

c)  The  synoptic  scale  departure  of  the  wind  from  geostrophic 
equilibrium  (the  sum  of  (a)  and  (b)). 

d)  The  smaller  scale  structure  of  the  wind  observed  by  the 
conventional  radiosonde. 

There  are,  as  well,  a  number  of  interesting  by-products  of  this  investigation. 
In  the  oourse  of  the  study,  statistics  on  the  magnitude  of  the  first  and  second 
horizontal  derivatives  of  the  field  of  isobaric  height  were  computed.  These 
results  should  be  useful  as  a  climtological  guide  in  constructing  or  testing 
models  of  large  scale  atmospheric  flow.  The  poleward  transport  of  moment  in 
was  found  to  occur  in  a  systematic  pattern  of  accelerations,  a  result  which 
has  such  broad  implications  for  models  of  the  general  circulation,  it  is 
discussed  in  a  separate  section  of  this  report, 

2.  Radionale  for  the  definition  of  synoptic  scale. 

The  spacing  of  upper  air  observing  stations  ultimately  will  set  a  limit  to 
the  horizontal  dimensions  of  systems  which  can  be  represented  on  the  weather 
map.  Although  the  spacing  of  upper  air  observing  stations  is  irregular,  the 
typical  spacing  in  continental  United  States  is  from  250  to  350  km.  The  present 
network  appears  capable  of  delineating  systems  with  horizontal  dimensions  in 


-  3  - 


excess  of  1000  km,  an  arbitrary  figure  which  is  often  used  as  the  limit 
of  synoptic  scale,  Ever,  an  isolated  jet  stream  which  has  a  transverse  di  Tension 
of  500  km  (e.g,  Endlich  and  McLean,  1965)  would  fall  within  this  definition  of 
synoptic  scale. 

Primarily  due  to  studies  nude  by  Buell  (1957,  1958,  1960),  it  is  known  that 
the  typical  dimensions  of  upper  air  systems  in  the  pressure  field  (the  pressure 
field  determined  from  the  radiosonde  observation  and  the  hydrostatic  computation 
of  pressure  height)  are  greater  than  2000  km  and  only  a  small  fraction  of  the 
variance  of  isobaric  height  can  be  associated  with  systems  of  smaller  scale. 
Studies  of  wind  structure  shew  that  the  variance  of  the  radiosonde  wind  must 
be  distributed  over  a  wider  range  of  scales  (Eddy,  1964)  although  the  spectral 
analysis  of  a  vector  quantity  requires  more  information  than  a  scaler  variable 
and  details  of  the  wind  structure  are  lacking. 

As  the  subsequent  analysis  will  show,  the  variance  of  isobaric  height 
is  90  nearly  restricted  to  larger  scale  systems,  the  analysis  of  the  pressure 
gradient  force  at  synoptic  scale  requires  little  smoothing  or  filtering;  the 
problem  is  really  one  of  interpolation  with  only  a  minimal  amount  of  smoothing. 
The  horizontal  contour  gradients  were  derived  by  a  process  of  numerical 
interpolation  and  smoothed  to  eliminate  systems  smaller  than  1000  km.  The 
pressure  gradient  force  determined  in  this  way  conveniently  defines  the 
synoptic  soale  pressure  systems.  The  wind  field  was  not  smoothed  on  this 
same  soale*  Instead,  a  ooeparisan  of  the  force  and  wind  field  was  accomplished 
by  interpolating  the  contour  gradients  at  the  location  of  a  radiosonde  wind 
observation.  The  scale  analysis  of  the  wind  field  was  then  nude  by  comparing 
the  dimensions  of  the  wind  system  to  the  dimtnsions  of  the  synoptic  soale 
pressure  field. 
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3.  The  nunerical  interpolation  of  the  pressure  gradient  force. 

In  oontrast  to  the  conventional  meteorological  analysis  of  the  field  of 
isobaric  height  which  is  used  to  interpolate  a  value  of  the  height  at  fixed 
grid  points,  the  purpose  of  this  analysis  is  to  interpolate  (as  it  finally 
de  sloped)  all  the  first  and  second  order  horizontal  derivatives,  five  derivatives 
in  all,  at  a  point  where  a  wind  observation  is  made.  The  irregularity  of  station 
spacing  as  well  as  observational  error  and  small  scale  variability  make  it 
necessary  to  interpolate  the  derivatives  from  a  larger  set  of  observations  than 
the  minimum  required  to  specify  five  derivatives.  After  some  experimentation 
with  interpolation  polynomials  of  various  degree  including  the  third,  a 
weighted  polynomial  of  second  degree  fitted  to  ten  observations  of  isobaric 
height  was  used  in  this  study.  Symbolically,  coefficients  of  the  polynomial 

hp  =  A2X^+A2xy  +  Ajy2'fA^X'*,Agy*Ag  (2) 

were  calculated  for  each  observing  period  by  minimizing  the  expression 

10  2 

£  thoi  -  (w<p  '  -J-  )  1 

i=l  R 

The  quantity  hQ^  is  the  isobaric  height  observed  at  the  ith  station  and  h^ 
the  height  at  that  station  predicted  by  the  interpolation  polynomial*  The 
origin  of  the  polynomial  wee  located  at  the  point  for  which  the  derivatives 
are  to  be  interpolated  and  the  exponent  of  the  error  type  weighting  function 
is  tire  square  of  the  distance  firm  the  observation  to  the  point  of  inter¬ 
polation  \  r^2  *  x.2  ♦  y.  2.  The  weighting  scale  parameter  R  was  set  to  equal 

1.  the  position  oFthe  observing  stations  were  calculated  in  a  Cartesian 
coordinate  system  with  the  axi  j  located  on  the  central  observing  station  of 
the  given  net.  So  that  the  analysis  could  be  readily  compared  with  conventional 
meteorological  charts,  the  c* xulation  was  performed  using  a  mathematical 
specification  of  a  Lambert  Conformal  Conic  projection  with  standard  parallels 
at  30°  and  6C°  latitude.  The  station  location  was  entered  to  the  nearest  minute 
of  latitude  and  longitude  according  to  the  Weather  Bureau  publication,  Key  to 
Meteorological  Records  Documentation  Mo.  5.21,  History  and  Catalogue  of  Uppe** 

Air  Data  for  the  Period  1946-1960,  (1964). 


500  km.  The  least  square  fit  with  this  weighting  factor  has  the  effect  of 
placing  primary  weight  on  the  four  or  five  observations  lying  with  a  radius 
of  500  km  and  of  also  smoothing  the  derivatives s  primarily  the  second 
derivatives,  an  a  scale  of  approximately  1000  km. 

The  second  degree  polynomial  has  frequently  been  used  in  objective  ’weather 
map  analysis,  most  recently  by  Corby  (1961).  It  is  worthy  of  reiteration, 
however,  that  in  almost  every  instance,  the  coefficient  A$  in  (2)  is  the 
value  interpolated.  In  the  present  analysis,  the  coefficients  A^  through 
are  the  end  product  of  the  interpolation  scheme. 

Coefficients  Aj^  and  A^  are*  of  course,  proportional  to  the  components  of 
the  horizontal  pressure  gradient  force.  If  the  coordinates  of  the  polynomial 
are  rotated  so  that  the  x'  axis  is  aligned  with  the  grostrophic  wind,  the 
first  three  coefficients  are  proportional  to  the  seoond  derivatives  and  will 
have,  as  well,  a  meteorological  interpretation  that  is  readily  recognized. 

The  coefficients  in  the  rotated  system  are  indicated  with  a  prime.  With 
the  appropriate  scaling  factors,  they  will  have  the  following  interpretation: 


2A'3 

a2h 

2  3y’2  * 

H  , 

V 

The  geostrophic  wind  shear, 
positive  for  cyclonic  shear. 

2  A'i 

a2h 

8  ax'2  * 

VgK„ 

The  geostrophic  wind  speed  x 
curvature  of  contours,  positive 
for  cyclonic  curvature. 

A'2 

a2h 

s  3xT  iy' 

vg*s.  • 

The  geostrophic  wind  speed  x 
curvature  of  orthogonals,  positive 
for  diverging  contours. 

4.  A  comparison  of  the  magnitudes  of  synoptic  geostrophic  winds  and  the 
observed  wind. 


Computations  were  made  for  seven  nets  of  ten  ot serving  stations  centered 
on  the  radiosonde  wind  observations  at  Flint,  Dayton,  Columbia,  Topeka,  Amarillo, 
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Shreveport  and  Montgomery,  tor  all  available  observations  in  the  winter  season 
December  1958  through  February  1959,  This  made  a  total  of  1200  comparisons  of 
geostrophic  wind  and  with  accompanying  wind  observation  at  the  500  mb  level 
and  1123  comparisons  at  the  3^r  mb  level  (of  a  possible  1274  data  periods). 

The  statistics  for  the  observed  wind  were  found  .o  be  reasonably  consistent 
with  published  values  based  cn  longer  periods  of  record  so  the  data  is  believed 
to  provide  a  representative  sample  for  the  winter  season.  While  some  variation 
in  the  statistics  was  observed  for  different  stations,  it  will  probably  be 
necessary  to  examine  data  far  several  seasons  to  establish  their  significance 
and  for  greater  generality  only  average^  statistics  for  the  five  stations 
with  best  records  are  presented  in  Table  I. 

Table  I  Averaged  Means  and  Variance  of  the  Components  of  the  Geostrophic 
Wind  and  the  ( jserved  Wind. 

5 CO  mb  300  mb 


Geostrophic 

Winds 

mean 

m  sec”^  1 

-variance  2 
o  |  nr  sec” 

mean  , 
m  sec”1 

-variance  2 
o  ,  tit  sec 

V  (Zonal) 
gx 

24.12 

93,2 

36,32 

201.6 

V  (Meridional) 

gy 

»3. 86 

126.5 

-4.48 

261.8 

Radiosonde  Winds 

V  (Zonal) 

X 

22,84 

112.8 

34.73 

275.7 

Vy  (Meridional) 

-3,10 

159.7 

-2.70 

342.7 

1c ,  Ihis  average 

r*p;*aents  the 

average  of  the  seasonal  values 

for  five  stations 

Shreveport  and  Montgomery  wind  records  were  less  regular  than  the  other  stations 
and  these  data  were  therefore  omitted  from  the  average.  The  500  mb  average  is 
based  on  878  observations  while  at  300  mb  there  are  829  observations  so  that 
these  values  should  be  fairly  representative  of  the  winter  season  in  mid-U.  S. 
locations. 


The  mean  observed  wind  is  very  nearly  equal  to  the  near,  geos  trophic  wind. 
Since  the  difference  of  the  nears  for  the  individual  station  were  all  of 
similar  algebraic  sign,  the  departure  of  the  mean  wind  from  the  near 
geostrophic  wind  is  believed  significant  and  representative  for  nid-U.S. 

The  500  mb  mean  wind  is  about  1  m  sec"1  sub-geos trophic  but  aligned  with 
the  contours,  the  300  mb  wind  is  also  slightly  sub-geostrophic  and  directed 
across  the  contours  toward  lower  heights  with  a  component  of  the  order  of 
1  m  sec"1. 

The  comparison  of  observed  wind  variance  and  geostrophic  wind  variance 
is  more  significant.  The  actual  wind  variance  of  each  component  is  larger 
than  the  geostrophic  wind  component  by  about  251  at  500  mb  while  there  is 
an  even  greater  relative  increase  at  330  mb  It  will  be  the  purpose  of  the 
remainder  of  the  analysis  to  allocate  tliis  greater  variance  of  the  radiosonde 
wind  among  three  categories:  (1)  Synoptic  scale  accelerations,  (2)  motions 
of  smaller  scale  and  (3)  error  in  wind  measurement. 


5.  The  model  for  synoptic  scale  wind  variance, 

A  straightforward  method  for  the  analysis  of  variance  according  to  scale 
is  to  compute  the  autocorrelation  function  and  its  Fourier  transform.  The 
spatial  autocorrelation  function  of  even  a  two  dimensional  vector  is  a  tensor 
that  requires  four  times  as  much  information  for  its  specification  as  for  a 
scalar  autocorrelation.  The  autocorrelation  function  for  the  i  and  j  velocity 
components,  for  example,  would  be: 


are  velocity  components  measured  at  the  ends  of  the  position  vector  r.  It  is 
unlikely  that  the  present  radiosonde  network  can  provide  satisfactory  statistics 
on  a  function  of  this  generality.  Certainly  the  2300  observations  used  in  this 


-  8  - 


study  would  not  begin  to  specify  this  function  so  that  a  direct  spectral 
analysis  of  the  wind  data  did  not  appeal  possible.  Instead,  the  statistics 
for  the  field  of  pressure  gradient  foroe  were  found  to  fit  a  model  of  a  Gaussian 
distribution  of  isobarie  height  proposed  by  Buell  (1960),  This  model  then 
provides  a  filter  for  the  separation  of  synoptic  scale  motions  from  motions 
of  lesser  scale.  The  essentials  of  Buell's  model  for  the  synoptic  scale 
variance  of  isobarie  will  therefore  be  reviewed  here. 

The  distribution  of  variance  of  isobarie  height  is  reasonably  isotropic 
as  shown  by  nearly  circular  isopleths  of  the  autocorrelation  coefficients  at 
a  number  of  locations  at  middle  latitudes  (Buell,  195 8 j  Bertoni  and  Lund, 

1963),  Buell  suggested  that  the  isopleths  were  quite  well  approximated  by 
the  Gaussian  distribution  function 

Ph  (r)  s  exp  (-r2/2L2)  (3) 

where  ph  is  the  autocorrelation  function  of  isobarie  height  and  L  is  a  horizontal 
scale  parameter  which  is  of  the  order  of  1000  km  for  the  mid-troposphere  and 
at  middle  latitudes.  This  function  not  only  provides  a  fairly  good  description 
of  the  autocorrelation  isopleths  at  moderate  distanoes,  r**  L,  but  it  also 
possesses  a  number  of  properties  that  make  it  suitable  for  a  synoptic  scale 
autocorrelation  function,  A  synoptic  scale  autocorrelation  function  of  isobarie 
height  must  have  a  second  derivative  (finite  curvature)  at  the  origin,  a 
requirement  imposed  by  the  necessity  that  the  geostrophic  wind  variance  must 
exis  for  synoptic  flow,  The  curvature  of  the  Gaussian  function  at  the  origin 
is  determined  by  the  parameter  L  and  to  the  extent  that  the  wind  is  geostrophic 
the  ratio  of  the  variance  of  observed  wind  to  the  variance  of  isobarie  height 
determines  the  value  of  L.  The  scale  parameter  L  determined  from  the  ratio  of 
wind  variance  to  isobarie  height  variance  is  apt  ''Tximately  the  same  as  that 
obtained  by  fitting  the  autocorrelation  function  at  greater  radius  (compare, 
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for  example,  values  of  L  deduced  by  Thompson,  1957,  and  Buell,  1958). 

Buell  further  deduced  if  the  variance  of  isobaric  height  is  isotropic  and 
homogeneous  the  autocorrelation  tensor  p  ,  ..  for  synoptic  geostrophic^  winds 

©  J 

naast  be  oapable  of  specification  by,  at  most,  two  correlation  functions,  the 
correlation  of  longitudinal  components  and  the  correlation  of  the  transverse 


components.  As  in  turbulence  theory,  we  shall  use  the  symbols  f  and  g  for 

s  g 

the  longitudinal  and  transverse  correlations,  the  subscript  g  indicating  the 
correlations  refer  to  the  geostrophic  wind. 


The  condition  that  the  divergence  of  the  geostrophic  wind  must  vanish  leads 
to  the  following  expressions  for  the  correlation  functions.  The  longitudinal 
correlation  function  will  have  the  same  form  as  the  autocorrelation  function 


for  isobaric  height, 


fg  =  exp  (-r2/2L2)  (4) 

and  the  correlation  of  the  transverse  components  will  be 

r2 

gg  =  (l-r2/L2)  exp  (-  ^2  >  (5) 

Inspection  of  the  trace  components  of  the  velocity  oorrealticn  tensor 
for  both  the  geostrophic  and  actual  winds,  pjj(r)  and  P22^  shcwe^  that 
they  did  have  a  shape  consistent  with  the  functions  (4)  and  (5)  and  would 
be  approximated  by  the  same  function  but  with  a  coefficient  lass  than  one. 

The  longitudinal  correlation,  for  example,  would  be  better  fit  by  the  function 


fg  *  Kg  exp  <-r2/2L2),  where  Kg  *  1. 


<f») 


Of  course  a  real  autocorrelation  function  will  approach  the  value  one  at  the 
origin.  Since  the  closest  spacing  between  observations  is  250  km,  an  auto¬ 
correlation  of  the  form  (6)  sinply  indicates  a  loss  in  correlation  (1-K)  in 


tha  first  ooqputation  interval.  Put  in  terms  of  the  spectral  distribution  of 

3,'  'the  geostrophic  wind  is  really  by  definition  a  quantity  defined  at  synoptic 
scale  so  that  the  adjective  synoptic  is  superfluous  used  with  geostrophic  wind 
as  a  description  of  the  horizontal  pressure  gradient  force. 


-  10  - 

variance,  it  would  be  interpreted  to  mean  that  a  fraction  of  the  variance  K 
must  be  associated  with  motions  of  synoptic  scale  while  the  fraction  (l-K)raust 
be  either  produoed  by  error  or  is  variance  associated  with  lesser  scales*  The 
dimensions  of  these  systems  of  smaller  scale  are  not  precisely  determined 
except  that  they  must  be  less  than  twice  the  spacing  of  the  distance  between 
observing  stations*  Our  hypothetical  autocorrelation  function  would  therefore 
assume  a  gap  in  the  spectrum  of  variance  between  synoptic  and  lesser  scales* 

The  model  autocorrelation  function  adopted  to  represent  the  autocorrelation 
of  both  geostrophic  and  actual  winds  was  of  the  form  suggested  by  (6).  The 
coefficient  Kg  was  assured  to  be  the  same  far  both  longitudinal  and  transverse 
correlations  so  that  where  (6)  is  the  exoression  for  the  autocorrelation  function 
for  the  longitudinal  components,  the  autocorrelation  function  for  the  transverse 
velocity  components  of  the  geostrophic  wind  would  be  given  by 

gg  =  Kg  (l-r2/!*2)  exp  (-r2/2L2)  (7) 

The  observed  wind  autocorrelations  were  fit  by  the  same  type  of  function  but 
with  different  coefficients.  For  the  autocorrelation  function  for  the 
longitudinal  components  of  the  observed  wind  the  expression  f^Ky  exp(-r2/2L2) 
was  used,  where  the  subscript  v  indicates  the  quantity  applied  to  the  obeerved 
wind. 

Trace  components  of  the  correlation  tensor  were  computed  far  seven  pairs  of 
stations  far  both  the  500  mb  and  300  mb  data  for  both  geostrophic  and  for  actual 
winds.  A  single  scale  parameter  L  and  four  values  of  the  parameter  K  were 
determined  from  the  56  data  points  by  a  least  square  fit  of  the  Gaussian  functions 
(6)  and  (7).  Four  values  of  K  arise  because  it  was  assumed  that  each  level  might 
have  a  characteristic  value  of  Kg  and  Ky  associated  with  the  geostrophic  and 
actual  wind.  The  quality  of  the  analytic  expressions  of  the  type  (6)  and  (7) 


x 
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as  approximations  to  the  autocorrelation  fucntion  is  illustrated  by  a  comparison 
of  the  average  function  for  the  two  orthogonal  components  with  the  computed 

correlations,  also  averages  of  two  components. 

2 

1/2  (f  +g)  =  KCl-  777  )  exp  C-r2/2L2)  (8) 

2L 

It  is  possible  then  to  coup  are  each  set  of  observations  with  a  single  curve 

for  the  function  assumed  for  the  autocorrelation  function.  It  must 

be  confessed  too  that  averaging  tends  to  suppress  irregularities 

in  Idle  component  correlations  representing  true  anisotropies  in  the  wind 

systems.  In  Figures  1  and  2  the  best  fitting  average  correlation  function 

given  by  (8)  is  shown  with  the  average  correlation  of  the  trace  component  of 

the  correlation  tensor.  If  the  distribution  of  isobaric  height  were  truly 

isotropic,  then  1/2  (P^  ^  +  Pg  22^  would  be  equal  to  the  average  of  the 

correlations  of  the  longitudinal  and  transverse  components  »l/2(f_  +  g  )  • 

6  5 

The  individual  geostrophic  wind  correlation  points  lie  above  the  correlation 
of  the  actual  wind  in  every  case.  The  relative  magnitude  of  the  intercept  Kv 
and  Kg  obtained  by  the  least  square  fit  of  (8)  are  clearly  established  even 
though  there  is  considerable  scatter  about  the  analytic  curve.  In  each  correlation 
set  the  actual  wind  correlation  is  about  ten  peroent  less  than  the  geostrophic 
wind. 

6.  The  variance  of  the  synoptic  wind  and  the  geostrophic  wind. 

The  scale  parameter  L  related  to  the  horizontal  dimension  that  was  obtained 
from  the  fit  shown  in  Figures  1  and  2  was,  L  =  1071  km,  a  value  much  the  same 
as  those  computed  by  Buell  (1958) .  The  parameters  Ky  and  Kg  are  (as  discussed 
above  and  also,  for  example,  by  Eddy,  1964)  interpreted  as  that  fraction  of 
the  total  variance  associated  with  the  modified  Gaussian  autocorrelation  function 
which  has  been  assured  to  represent  the  variance  at  synoptic  scale.  The  equivalent 
spectrin  of  the  Gaussian  autocorrelation  function  is  well  known  and  drops  off 
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rapidly  in  the  wavelength  region  less  than  2L.  A  wavelength  parameter  of 
1000  km  implies  that  synoptic  scale  systems  have  dimensions  for  the  most 
part  greater  than  2000  km  and  systems  of  dimension  L  have  been  completely 
filtered  from  the  synoptic  variance  by  this  autocorrelation  function. 

The  values  of  the  parameters  Ky  and  Kg  for  the  500  and  300  mb  data 
obtained  by  a  least  square  fit  of  the  functions  of  the  type  given  by  (7) 
are  shown  in  Table  II,  The  resulting  division  of  variance  at  synoptic  sca^e 
and  at  lesser  scales  implied  by  the  analysis  is  also  shown.  In  all  cases, 
the  variance  in  Table  II  represents  the  variance  of  one  component  obtained  by 
averaging  ?onal  and  meridional  components. 

Table  II  Estimate  of  the  magnitude  of  synoptic  wind  variance  and  the 

variance  at  smaller  scale,  (synoptic  scale  dimension  parameter, 

L  =  1071  km). 


Fraction  of  total 
variance  in  synoptic 
scale  motion,  K 

Synoptic  scale 
variance  of  one 
component 
sec“^ 

Smaller  scale 
Variance  or  error 
Cl-K)  o2 
jit  sec 

500  mb  300  mb 

500  mb 

300  mb 

500  mb 

300  mb 

Geostrophic  wind 

0.940  0,932 

103.4 

215.9 

6.6 

15.8 

Radiosonde  winds 

0.817  0.821 

111.4 

254.0 

11, 3 

55.2 

The  rather  close  correspondence  of  the  values  of  Kg  and  Ky  obtained  from 
the  500  mb  and  the  300  mb  data  while  not  really  adding  validity  to  their 
interpretation  as  the  fraction  of  variance  in  synoptic  scale  has  been  of  some 
oorafort  to  the  analyst.  The  total  variance  of  the  observed  wind  is  so  much 
greater  than  the  total  variance  of  the  geostrophic  wind  that  despite  its  smaller 
fraction  of  synoptic  variance  the  variance  of  the  observed  wind  is  greater  than 
the  variance  of  the  geostrophic  wind  in  both  categories  of  scale. 
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The  geostrophic  wind  has  really  been  computed  at  synoptic  scale.  The 
estimate  of  variance  at  smaller  scale  or  error  far  the  geostrophic  wind  is  almost 
entirely  variance  produced  by  error  in  the  isobaric  height.  Indeed,  the  values 
(Table  II)  are  consistent  with  the  error  estimates  of  the  geostrophic  wind  made 
by  Endlich  and  Clark  (1963).  The  greater  variance  of  the  observed  wind  at 
synoptic  scale  as  compared  to  the  geostrophic  wind  is  assumed  to  be  associated 
with  the  geostrophic  departures  at  synoptic  scale  or  the  synoptic  scale 
accelerations,  in  Eq,  (1),  The  excess  varianoe  of  one  component  is 
8.0  m^sec”^  at  500  nib  and  28  m^sec“^  at  300  mb.  The  hypothesis  that  this 
difference  represents  the  geostrophic  departure  of  the  wind  at  synoptic  scale 
was  investigated  further  by  an  analysis  of  the  observed  geostrophic  departures. 

6.  The  synoptic  scale  departure  of  the  wind  from  geostrophic  equilibrium. 

An  independent  estimate  of  the  synoptic  accelerations  (the  accelerations 
measured  relative  to  the  coordinate  system  fixed  on  rotating  earth)  is  obtained 
by  an  analysis  of  the  departure  of  the  observed  wind  from  geostrophic  balance. 
Ideally  perhaps,  the  winds  would  be  smoothed  on  the  same  scale  as  the  pressure 
field  and  synoptic  accelerations  computed  from  the  difference.  In  this  study 
the  departures  of  the  observed  wind  from  the  synoptic  geostrophic  wind  were 
computed  and  the  departures  subsequently  analyzed  for  evidence  of  a  synoptic  scale 
variation. 

Despite  very  little  quantitative  evidence,  it  is  almost  a  meteorological  axiom 
that  the  geostrophic  departures  are  greater  in  the  direction  of  the  oon tours, 
the  longitudinal  direction,  thm  in  the  direction  across  or  transverse  to 
the  contour*.  Geostrophic  departures  were,  therefore,  computed  in  natural 
coordinates  es  well  as  in  zonal  and  meridional  corponents .  The  longitudinal  geo¬ 
strophic  departure  is  indicated  by  the  symbol  DT  and  oamputed  thus  DT*V  cos  3~Vg» 
where  ®  is  the  angle  made  by  the  wind  vector  with  the  oon  tours ,  positive 
toward  high  pressure.  The  transverse  geostrophic  departure  is  given 


1 


-  14  - 


by  Dn  s  V  sin  a  where  V  is  the  observed  wind  speed.  The  zonal  and  meridional 

components  are  simply  Dv  =  Vv-V_,  and  D  =  V  -V  positive  axis  toward  the 

x  X  gx  y  y  gy 

east  and  north,  respectively .  Average  values  for  the  geos  trophic  departures 
for  the  sane  five  stations  as  Table  I  (based  on  more  than  800  cases)  are 
given  in  Table  III. 


Table  III  Statistical  characteristics  of  the  departure  of  the  rawin sonde 
wind  from  the  synoptic  geos  trophic  wind. 


500  mb 

300 

mb 

Mean 
m  sec' 

o  m  sec 

Mean  . 
m  sec"1 

Variance 
o  m  sec 

Longitudinal  Component,  D 

0.97 

31.8 

1.16 

91.0 

Transverse  Component, 

-0.02 

18.4 

-1.06 

50.9 

Zonal  Component,  Dx 

1.28 

31.0 

1.70 

85.4 

Meridional  Conponent,  Dy 

-0.75 

18.3 

-2.44 

50.0 

The  magnitude  of  the  seasonal  mean  geoB  trophic  departures  of  the  order  of 
a  meter  sec**1  was  already  implied  when  the  mean  geoa trophic  wind  and  mean 
rawinsorvde  wind  were  entered  in  Table  I.  The  variance  of  the  longitudinal 
geos  trophic  departure  is  significantly  larger  than  the  variance  of  the 
transverse  component.  The  magnitude  of  these  computed  geostrophic  departures 
are  perhaps  best  put  in  perspective  by  the  diagram  in  Figure  3.  The  mean 
geostrophic  wind  vector  is  shown  for  the  500  mb  and  the  300  mb  levels  and 
the  ellipse  representing  one  standard  deviation  of  the  departures  of  the  wind 
from  the  geostrophic  wind  is  drawn  about  the  terminus  of  the  mean  geostrophic 
wind  vector.  Typical  cross  contour  angles  of  the  wind  (a)  are  indicated  by 
the  vector  drawn  to  the  boundary  of  the  ellipse.  Thus  even  though  the  standard 
deviation  of  the  zonal  component  of  the  wind  is  almost  20  knots  at  300  mb,  the 
typical  cross  contour  angle  of  10°  is  within  synoptic  experience  since  the  winds 
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are  only  reported  to  the  nearest  10°  in  direction4. 

Since  one  would  like  to  attach  physical  significance  to  the  observation 
that  the  variance  of  the  longitudinal  geos  trophic  departure  is  50%  to  75% 
lazier  than  the  variance  of  the  transverse  component,  it  would  aid  this 
interpretation  if  the  variance  in  the  zonal  and  meridional  direction  were 
more  nearly  equal.  The  zonal  departure,  however,  is  the  same  magnitude  as 
that  of  the  longitudinal  component  and  the  meridional  component  is  'die  same 
magnitude  as  the  transverse  component.  Apparently  because  the  wind  is  so 
nearly  zonal,  it  is  inpossible  with  these  data  to  resolve  the  precise 
orientation  of  the  ellipse  shown  in  Figure  3,  One  is  not  able  to  distinguish 
between  the  possibility  that  the  observed  difference  is  primarily  a  zonal- 
meridional  difference  or  whether  it  is  better  related  to  the  direction  of 
the  synoptic  flow. 

To  investigate  whether,  in  fact,  the  geostrophic  departures  possessed 
any  component  at  synoptic  scale  the  autocorrelation  of  geostrophic  departures 
were  computed  for  the  same  seven  pairs  of  stations.  The  average  correlation  of 
the  two  components  of  geostrophic  departure  are  shown  as  a  function  of 
separation  in  Figure  4,  The  300  mb  data  is  sufficiently  regular  so  that  a 
curve  of  the  form  of  (5)  oan  readily  be  fitted  to  the  data.  The  curve 
shown  has  the  intercept  K  *  0,52,  and  a  horizontal  dimension  parameter  L  = 

960  km,  very  similar  to  that  found  for  the  geostrophic  wind. 

The  500  mb  geostrophic  departure  correlations  are  more  irregular  and 

do  not  clearly  exhibit  the  same  horizontal  scale  as  the  300  mb  departures. 

It  may  be  appropriate  to  consider  the  extrapolated  intercept  an  upper  estiirute 

of  the  fraction  of  geostrophic  departures  at  synoptic  scale.  The  estirates  of 

4.  Tn  the  climatological  data  used  in  this  study  the  wind  directions  were 
reported  to  1®,  It  may  be  of  some  interest  that  two  large  geostrophic 
departures  were  ultimately  traced  to  errors  remaining  in  the  checked 
climatological  wind  data. 
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synoptic  geos  trophic  departures  obtained  from  the  intercept  and  the  variance 
of  the  geos  trophic  departures  (K  o  )  for  one  component  are  as  follows: 

At  300  mb  36.  m^sec  ^ 

At  500  mb  9.5  m^sec"^ 

These  are  values  slightly  larger  but  consistent  with  the  estimate  obtained 
fnan  the  excess  synoptic  variance  of  the  rewinscnde  winds.  Geos  trophic 
departures  are  larger  at  300  mb  than  at  500  mb  and  a  larger  fraction  of  the 
variance  of  the  geos  trophic  departure  has  the  dimension  of  synoptic  scale. 
Stating  the  result  in  a  slightly  different  way,  the  synoptic  scale  wind  at 
500  mb  is  more  nearly  geos  trophic  and  has  as  a  consequence  less  divergence  than 
the  wind  at  300  mb. 

7.  Patterns  of  synoptic  accelerations. 

The  observation  that  even  with  highly  curved  streamlines  and  rapidly 
changing  pressure  patterns  the  wind  vector  remains  aligned  with  the  contours 
has  led  to  many  proposals  for  direct  relations  between  the  accelerations  and 
the  pressure  field  or  what  is  equivalent,  proposals  for  mechanisms  of  adjustment 
of  the  wind  to  geos  trophic  balance.  The  equations  of  motion,  on  the  other  hand, 
predict  that  accelerations  must  depend  on  the  initial  conditions  and  a  pattern 
of  accelerations  related  to  the  p^-tem  of  pressure  gradient  force  may  be 
equivalent  to  requiring  a  synoptic  frictional  force  of  essentially  the  same 
magnitude  as  the  synoptic  accelerations.  The  above  estimates  of  the  magnitude 
of  synoptic  accelerations  (that  a  third  to  a  half  ot  the  observed  wind  departure 
from  ge oe trophic  may  be  of  synoptic  scale)  enoouraged  the  test  of  r-crae  features 
of  the  models  of  synoptic  acceleration. 

Althoufh  the  derivation  of  models  of  acceleration  differ  (and  in  some  cases 
do  yield  conflicting  predictions)  there  is  a  fundamental  uniformity  in  the 
prediction  that  the  first  approximation  to  the  synoptic  acceleration  will  be 
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given  by  the  acceleration  of  the  geostrophic  wind,  dV/dt  3  dVg/dt  (e.g.  Chamey, 
1955;  Bykov,  1962;  Petterssen,  1956;  Amascm  et  al,  1262),  The  gradient  wind 
approximation  can  also  be  construed  to  be  of  this  form.  The  only  terms  of 


that  can  be  evaluated  on  an  isobanc  chart  is  the  advective  component 

V  *  Vp  and  these  terms  were  computed  and  compared  with  the  observed 

geostrophic  departures.  The  geostrophic  departures  DT  and  Dn  were  correlated 

with  the  components  of  the  geostrophic  wind  acceleration  [(V  cos  a)  (V^  ] 

and  {(V  oos  a)  (V  K  )}  and  the  results  given  in  Table  IV.  The  relation 

g  n 

dV  dVa 

U|r  =  -K&.  implies,  among  other  things,  that  the  wind  should  be  sub- geostrophic 
with  cyclonically  curved  contours  and  be  directed  toward  high  pressure  with 
converging  contours. 


Table  IV  The  correlation  of  the  departure  of  the  wind  from  geostrophic  with 
acceleration  of  the  geostrophic  wind. 


590  mb 

300  mb 

dt 

VS. 

[<V 

oos 

a) 

<V8 

VJ 

9.061 

0.316 

vs. 

[(V 

cos 

a) 

tV8 

V1 

0.138 

0.420 

Dt 

vs. 

[(V 

cos 

3  ) 

% 

0.071 

-3.006 

Dn 

vs. 

[(V 

oos 

a) 

K,)J 

0.076 

0.052 

The  first  two  correlations  are  those  anticipated  by  the  relation 
dV  dV„ 

=  -jJk  ;  the  second  set  were  computed  simply  as  a  measure  of  the  noise 
level  in  these  data. 

Neither  of  the  correlations  with  geostrophic  departure  with  the  terms  of 

dV„ 

-Ja  are  very  significant  at  the  590  mb  level  but  at  30G  mb  both  relations  are 
almost  as  high  as  cool-  be  expected  in  view  of  the  previous  estimates  of  snail 
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scale  variability  of  the  geostrophic  departures c  Despite  the  fact  that  is 

larger  than  Dn,  the  correlation  of  Bt  with  the  contour  curvature  term  is 

consistently  less  than  the  correlation  of  Dn  with  the  term  involving  the 

convergence  of  the  contours.  The  slight  difference  in  the  correlations  does 

not  properly  demonstrate  the  rather  marked  difference  in  the  quality  of 

the  two  relations.  The  two  sets  of  data  for  300  mb  for  Dayton  are  plotted 

on  scatter  diagrams  to  illustrate  this  point.  In  Figure  5,  the  transverse 

geostrophic  departures  are  plotted  against  the  gecstrophic  wind  acceleration 

associated  with  contour  convergence.  The  line  of  best  fit  very  nearly 

,j\j  dV„ 

corresponds  to  a  line  ^  ,  In  Figure  6,  which  is  a  scatter  diagram 

for  the  relation  of  the  longitudinal  geostrophic  departure  plotted  versus  the 

geostrophic  wind  acceleration  associated  with  contour  curvature,  it  is  clear 

the  relatively  high  correlation  has  not  been  achieved  by  a  normal  distribution 
^  , .  dV 

aoout  the  line  , 

8.  The  magnitude  of  the  synoptic  friction  tern. 

That  the  synoptic  accelerations  can  be  approximated  by  the  acceleration 

dV 

of  the  geostrophic  wind  has  been  used  to  justify  the  basic 

assumption  employed  in  the  balanced  wind,  (V  *V)  =  0,  The  balanced 

wind  has  also  been  termed  a  "frictionless”  approximation  to  the  flow  +  -*■ 

d  V 

CPanofsky,  1957),  However,  it  is  readily  seen  that  the  observation  eft  3t 

may  with  equal  validity  be  interpreted  to  mean  that  the  frictional  force 
must  be  the  same  magnitude  as  the  synoptic  accelerations.  If  in  Cl)  the 
geostrophic  wind  relation  is  substituted  for  the  horizontal  pressure  gradient 
force,  cne  obtains 


§  =  f  X  x  cvB  -  V)  +  r„  . 
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Differentiation  of  this  expression  with  respect  to  time  yields 


is  the  relation  suggesting 


If  the  vector  standard  deviation  of  the  synoptic  acceleration  is  used 
as  a  measure  of  the  magnitude  of  the  acceleration  and  the  frictional  force, 
the  relative  magnitude  can  best  be  illustrated  as  a  fraction  of  the  average 
pressure  gradient  farce  at  that  level.  At  500  mb  this  ratio  will  be 
o  dV/dt  (g  V  p  h)~*  =  0*16  and  at  300  mb,o  dV/dt  (g  v  p  h)“*  =  0.22. 

It  should  not  be  inferred  from  this  that  the  frictional  force  is  directly 
proportional  to  the  pressure  gradient  force  for  it  is  clearly  asymmetricaJL 
in  just  the  same  way  as  the  geos  trophic  departures  are  asymetrical.  The 
meridional  or  transverse  frictional  foroe  is  perhaps  50%  greater  than  the 
zonal  or  longitudinal  frictional  force.  For  the  same  seasons  discussed 
earlier,  it  is  not  possible  to  distinguish  between  the  merits  of  a  model  of 
the  frictional  foroe  proportional  to  the  geostrophic  departure  (assumed,  fo^ 
example,  by  Kao,  1963)  or  one  with  a  zonal  meridional  asymmetry  as  assumed 
by  Davies  and  Oakes  ,1962,  either  form  vould  be  consistent  with  these  data. 

The  scatter  diagram  of  Figure  6  is  cited  as  evidence  for  a  transverse 
frictional  force  of  the  order  of  20%  of  the  mean  pressure  gradient  foroe  at 
300  mb.  Although  the  correlation  between  DT  and  dV/dt  is  significant,  the 
data  points  spread  out  in  a  V-shaped  region  so  that  the  line  dVg/dt  is  only 
a  rough  bound  to  the  sub-geostrophic  departures  but  large  super-gecstrophic 
winds  are  also  found  with  cyclonically  curved  contours.  Since  the  transverse 
geostrophic  departures  are  snail  and  independent  of  the  longitudinal  geostrophic 
departure,  it  seems  clear  that  the  alignment  of  the  wind  with  the  contours  can 
only  be  maintained  by  invoking  forces  other  than  the  synoptic  pressure  gradient 
foroe. 
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9.  The  poleward  transport  of  angular  moment  un. 

It  is  well  known  that  the  sink  of  angular  momentum  at  the  surface  of 
the  earth  in  raid- latitudes  must  be  fed  aloft  by  a  transfer  of  momentum  from 
lower  latitudes.  The  transport  of  angular  momentum  computed  from  the 
correlation  of  the  zonal  and  meridional  wind  components  is  quite  variable  and 
dependent  upon  station  location  but  at  mid- continental  locations  in  the 
central  states,  it  would  be  expected  to  be  positive  at  both  500  and  300  mb 
levels.  The  correlation  of  the  observed  wind  components  was  in  fact  positive 
at  all  stations  and  the  winter  five  station  average  correlation  coefficient 
vx  Vy  ovx_1  ovy”1,  was  computed  to  be  0,081  at  500  mb  and  0,307  at  300  mb. 
These  figures  are  not  necessarily  representative  of  the  transport  angular 
momentum  at  these  latitudes  but  are  believed  representative  of  this  region 
for  the  winter  season  under  investigation.  The  correlation  coefficient  for 
the  geostrophic  wind  components  was  also  computed  and  while  the  absolute 
magnitude  of  the  transport  may  be  in  doubt,  it  was  found  that  the  northward 
transport  of  angular  momentum  by  the  observed  wind  exceeded  the  northward 
transport  by  the  geostrophic  wind  at  all  stations  and  at  all  levels.  The 
corresponding  five  station  average  correlation  coefficients  for  the  geostrophic 
wind  were  -0.031  and  0,104,  It  may  be  concluded  therefore  that  for  this 
region,  the  poleward  transport  of  momentum  occurred  largely  in  the  ageo- 
s trophic  component  of  the  wind. 

The  transport  of  angular  momentum  was  investigated  further  at  the  300  mb 
level  where  the  difference  between  the  correlation  coefficients  appeared  to 
be  most  significant.  Simply  by  dividing  the  cases  into  either  northerly  or 
southerly  flow,  it  was  found  that  on  the  average  north  winds  were  almost 
parallel  to  the  contours  and  slightly  sub-geostrophic  but  most  of  the 
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difference  in  the  correlation  was  produoed  by  only  20%  of  the  cases  of 

southerly  flow.  In  southerly  flow  there  is  an  average  geostrophic  departure 

at  all  stations  j  the  wind  speed  is  super-geostrophic  by  several  meters  per 

second  and  directed  across  the  contours  toward  lower  isobaric  height  at  an 

angle  cf  almost  10°,  Since  southerly  flow  at  300  mb  is  usually  associated 

with  upward  motion  and  positive  wind  shear  these  observed  geostrophic 

dV  dVg 

departures  are  consistent  with  the  relation  jjr  =  ,  the  geostrophic 

acceleration  in  this  case  being  the  vertical  advection  of  the  geostrophic 
wind. 

In  summary,  the  northward  transport  of  angular  momentum  at  300  mb  is 
largely  produced  by  geostrophic  winds  in  southerly  flow.  The  wind  is  on  the 
average  slightly  super-geostrophic  and  directed  toward  lower  isobaric  height. 

10,  Summary  of  results. 

The  magnitudes  of  the  velocity  fluctuations  and  accelerations  determined 
in  this  study  are  summarized  by  the  estimate  of  the  standard  deviation  of  the 
velocity.  Standard  deviations  of  various  quantities  are  given  in  equivalent 
velocity  fluctuations  with  units  of  meters  sec“^  in  Table  V,  The  mean  wind 
speed  at  each  level  is  tabulated  far  reference.  In  view  of  the  sample  size 
the  figures  have  been  rounded  to  the  nearest  whole  number. 


Table  V  Estimated  standard  deviation  of  the  velocity  in  meters  sec"^-, 
mid-Uc  S.,  winter  season - 

500  mb  300  mb 

Mean  Wind  Speed  : j - J  34,8 

Standard  deviation 

Observed  geostrophic  departures  all  scales 

Dt  longitudinal  component  6.  10, 


transverse  component  4  7 

Synoptic  scale  acceleration  of  one 

component,  equivalent  velocity  fluctuation,,  3  6 

Small  scale  ve  loci  tv  fluctuation  of 

radiosonde  wind  (includes  error)  3  7 

Error  of  geostrophic  wind  2  4 


The  small  scale  fluctuation  of  velocity  and  the  synoptic  acceleration 
are  about  the  same  order.  The  synoptic  scale  divergence  will  therefore  also 
be  about  half  the  observed  divergence  a  result  consistent  with  a  study  by 
Eddy  (1964),  The  synoptic  scale  vector  acceleration  is  almost  20%  as  large 
as  the  horizontal  pressure  gradient  force,  This  figure  of  20%  must  also 
represent  an  upper  ragnitude  for  the  frictional  force  in  the  free  atmosphere. 

Appendix:  Seme  characteristics  of  the  synoptic  pressure  field. 

Our  ability  to  filter  synoptic  scale  systems  from  the  background  of  smaller 
scale  effects  and  error  is  largely  dependent  cn  the  fit  of  a  single  parameter 
Gaussian  function  to  the  autocorrelation  function  for  iso  baric  height.  The 
aurves  of  Figures  1  and  2  illustrate  the  fit  of  this  function  at  large  distance 
from  the  origin  and  show  that  the  scale  parameter  L  =  1071  km  is  a  correct  order 
of  magnitude  fit  at  both  the  330  mb  and  500  mb  levels,  Buell  (1954)  noted 
that  if  the  horizontal  gradient  of  the  variance  of  isobaric  height  is  srrall 
(a  condition  required  for  homogeneity)  the:  the  curvature  of  the  correlation 
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function  of  isobaric  height  at  the  origin  will  be  given  by  the  ratio  of  the 
variance  of  the  geostrophic  wind  to  the  variance  of  isobaric  height.  In  the 
notation  of  the  text 


3  ph(0) 


3x 


2 

3  Ph(0) 


=  -Cf/g)2 


o2  Vgx  +  a  Vgy  -j 

ah2 


(la) 


If  p(r)  =  exp  ( 
parameter  L, 

Values  of  L  obtained  from  the  relation  (2a)  are  consistent  with  the  value  of 
1071  km  obtained  in  the  least  square  fit  of  data  at  large  distances  from  the 
origin. 

Sane  anisotropies  in  the  field  of  isobaric  height  were  noted,  however, 
as  systematic  differences  in  the  magnitude  of  the  second  derivatives.  At 
300  mb  in  particular  the  anisotropy  could  be  related  to  the  position  of  the 
jet  stream  and  some  of  the  variation  in  the  statistics  with  station  are 
therefore  considered  significant.  Statistics  on  the  seocnd  derivatives  of 
isobaric  height  at  the  five  stations  with  most  complete  reoords  are  given 
in  Table  (IA).  The  mean  value  of  the  seocnd  derivatives  for  the  winter 
season  are  tabulated  with  the  standard  deviation  in  mits  of  equivalent 
vorticity  or  shear.  » 


2L 


),  then  the  curvature  at  the  origin  determines  the  scale 


Ls  2l/2B/f  «h  V  *  0  vgy  * 
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The  seasoned  average  of  the  second  derivatives  are  small  at  all  stations. 
A  small  cyclonic  curvature  and  a  slight  confluence  of  the  mean  contours  is 
observed  at  all  stations  as  indicated  by  t)ie  positive  values  of  the  mean 


derivative 


32h 


— k  and  the  negative  values  of  the  mean  derivative 

3v  * 


32h 


<*xr  3x3y' 

This  pattern  of  oon tours  also  appears  in  published  mean  seasonal  contour  charts 

at  these  locations* 

The  distribution  of  values  of  t^e  derivatives  measured  in  natural 
coordinates  closely  approaches  a  normal  distribution  so  that  the  standard 
deviation  provides  a  good  estimate  of  the  representative  magnitudes.  For 
example,  in  only  one  case  did  the  anti -cyclonic  geos  trophic  shear  exoeed 
the  value  of  the  Coriolis  parameter, a  value  approximately  three  standard 
deviations  from  the  mean*  The  principal  anisotropy  shown  by  these  data 
is  the  tendency  towards  a  zonal  or  longitudinal  striation  as  shown  by  the 
consistently  greater  magnitude  of  the  geostrophic  shear  as  compared  to  the 
longitudinal  second  derivative  associated  with  oon  tour  curvature.  Only  at 

Amarillo  is  there  approximate  isotropy  in  the  seoond  derivatives. 

32h 

The  smeller  magnitude  of  the  cross  derivative,  jxTyyi-"  ,  does  not  indicate 
anisotropy.  It  is  instead  a  measure  of  the  ratio  of  the  amplitude  of  synoptic 
scale  velocities  to  the  speed  of  the  zonal  current.  In  any  case,  these 
derivatives  have  almost  the  same  magnitude  at  a  given  level  at  all  stations* 

The  vorticity  varianoe  is  approximately  that  given  by  the  sum  of  the  varianoe 
of  each  derivative  indicating  almost  no  correlation  in  the  component  seoond 
derivatives.  The  vorticity  increases  by  at  least  fifty  percent  from  500  to 
the  300  mb  level  and  is  remarkably  uniform  with  station  location. 

The  magnitude  of  the  seoond  derivatives  can  be  illustrated  in  another  fashion 
by  a  computation  of  typical  values  of  the  oontour  radius  of  curvature  and  radius 
of  curvature  of  the  orthogonals.  A  radius  of  curvature  of  the  contours  was 
ocnputed  from  the  ratio  of  the  mean  geostrophic  wind  speed  to  the  standard  deviation 
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of  the  longitudinal  derivative 


This  is  in  effect  another  measure  of  the  size  of  pressure  systems;  values 
of  the  radius  of  curvature  of  the  oontours  calculated  at  these  stations  is 
given  in  Table  I  LA.  The  values  vary  with  the  station  location  but  are 
almost  constant  with  altitude. 

Values  of  the  orthogonal  curvature  showing  convergence  or  divergence  of 
the  oontours  are  computed  in  similar  fashion  from  the  standard  deviation  of 
the  cross  derivative.  The  orthogonal  curvature  is  almost  constant  with 
altitude.  Typical  values  of  orthogonal  curvature  are  given  in  Table  IIIA, 

In  conclusion  the  characteristic  of  the  synoptic  scale  field  of  isobaric 
height  are  fairly  well  summarized  as  isotropic  disturbances  on  a  zonal  current 
with  a  horizontal  scale  greater  than  10CO  km.  The  scale  of  the  systems  is 
not  greatly  different  at  S00  and  300  mb.  The  major  anisotropy  is  the 
longitudinal  stziaticn  associated  with  what  might  be  called  the  jet  stream. 
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Table  IIA  Typical  values  of  radius  of  curvature  of  the  contours,  km. 
Another  measure  of  the  size  of  pressure  systems. 


SOO  mb 

300  mb 

Flint 

1540  km 

1500  km 

Dayton 

1600 

1680 

Columbia 

1240 

1300 

Topeka 

1190 

1090 

Amarillo 

950 

950 

Table  IIIA  Typical  values  of  the  radius  of  curvature  of  oontour 

orthogonals,  km.  A  measure  of  the  degree  of  confluence 
of  the  contours. 


500  mb 

300  mb 

Flint 

2500  km 

2560  km 

Dayton 

2290 

2710 

Columbia 

1930 

1930 

Topeka 

1990 

1970 

Anariilo 

2240 

1920 
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FIGURE  CAPTIONS 


Fig.  1  Average  autocorrelation  of  traoe  components  of  the  geos  trophic 
wind  (solid  circles)  and  rawinsonde  wind  (open  circles)  for 
the  50C  mb  level,  Jhe  curve  is  the  least  square  fit  of  the 
function,  K(l-r  /2L* )  exp  (  -  r2  ). 

2L2 

Fig.  2  Average  autocorrelation  of  trace  components  of  the  geostrophic 
wind  (solid  circles)  and  rawinsonde  wind  (open  circles)  for 
the  300  mb  level ,  The  curve  is  the  least  square  fit  of  the 
function,  K(l-r*/2L2)exp  (  _  ). 

2L2 

Fig.  3  Illustration  of  representative  total  departure  of  the  wind  from 
geostrophic  computed  in  this  study.  The  mean  geostrophic  wind 
vector'  is  plotted  to  scale  and  ellipses  of  semi  major  and  semi 
minor  axes  of  one  standard  deviation  of  the  zonal  and  meridional 
geostrophic  departures  are  drawn  about  the  terminus  of  each 
vector.  The  cross  contour  angle  should  lie  within  the  angle  a 
about  two  thirds  of  the  time. 


Fig.  4 


Fig.  5 


Average  autocorrelation  of  orthogonal  components  of  the  departure 
of  the  wind  from  geostrophic.  Solid  circles  300  mb  data,  open 
circles  500  mb  data.  The  curve  is  the  fit  of  the  function 


to  the  300  mb  data. 


Scatter  diagram  of  the  relation  between  the  transverse 
geostrophic  departure  and  the  longitudinal  component  of 
geostrophic  acceleration.  Example  of  Dayton  data  300  mb. 
The  broken  line  represents  the  relation  predicted  by  the 
assumption  -* 
d  V 

Ht  = 


Fig.  6  Scatter  diagram  of  the  relation  between  the  longitudinal 

geostrophic  departure  and  transverse  component  of  geostrophic 
acceleration.  Example  of  Dayton  data  300  mb.  The  broken 
line  represents  the  relation  predicted  by  tie  assumption 

dV 
(Jt 
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The  Vertical  Distribution  of  Dust 
by 

James  M.  Rosen 

Six  flights  have  been  rrade  with  a  coincidence  type  dust  counter 
similar  to  the  one  previously  described  (Ann.  Prog,  Report,  1964),  The 
results  of  two  of  these  flights  have  already  been  reported,  (Rosen,  1965, 

1966), 

The  first  flight  made  with  the  coincidence  dust  counter  is  a  good 
illustration  of  the  layered  character  of  the  stratosphere.  The  data 
are  shown  in  Figures  1  and  2.  The  region  between  the  tropopause  and  the 
ozone  maximum  is  frequently  characterized  by  dust  and  ozone  rivers  described 
by  Kroening  and  Ney  (1962),  On  this  sounding  there  appears  to  be  a  good 
correspondence  between  the  two  lower  dust  and  ozone  peaks.  At  19.5  km 
however,  there  is  an  ozone  river  but  no  corresponding  dust  river.  This 
observation  may  seem  in  conflict  with  the  hypothesis  that  the  dust  and 
ozone  are  transported  together  but  can  be  understood  by  referring  to  the 
lines  of  constant  mixing  ratio  on  the  sounding.  Since  the  mixing  ratio  of 
a  trace  constituent  is  preserved  in  vertical  displacements  a  river  will  not 
appear  in  a  sounding  for  that  constituent  in  a  region  where  the  mixing  ratio 
is  constant  with  altitude.  The  river  in  this  case  is  in  a  location  where 
the  dust  follows  a  oonstant  mixing  ratio  with  altitude  but  tnc  ozone  does 
not.  Thus,  only  the  o^one  profile  should  show  the  river.  Above  the  ozone 
maximum  where  both  the  dust  and  ozone  mixing  ratio  is  constant  with  altitude 
a  river  oannot  be  observed  in  a  sounding  for  these  constituents. 

A  casual  inspection  of  the  du_u  and  ozone  layers  at  about  17  km  would 
suggest  that  it  is  also  a  river,  but  upon  closer  examination  a  number  of 
peculiarities  appear  nonrally  not  associated  with  rivers.  First,  the  dust 
\M\mr  in  about  1  km  lower  than  the  ozone  layer  (the  resolution  is  about  1/6  km). 
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Second,  the  mixing  ratio  of  dust  in  this  layer  is  greater  than  at  any  point 
in  the  atmosphere  above  it.  In  addition,  the  size  distribution  in  this 
layer  is  considerably  different  from  that  of  the  surrounding  atmosphere 
(Figure  2). 

This  layer  is  perhaps  associated  with  the  radioactive  debris  from  the 
underground  explosion  in  Nevada  two  weeks  before  the  flight  or  from  the 
Chinese  explosion  six  months  before  the  flight, or  from  high  altitude  testing 
two  years  before  the  flight.  The  first  two  sources  seem  rather  unlikely  in 
view  of  the  fact  that  there  was  very  little  fallout  connected  with  them 
(Kuroda  et  al,  1965).  The  third  source  also  does  not  seem  very  likely 
because  there  is  over  two  years  between  injection  and  measurement.  However, 
thin  radioactive  layers  in  the  stratosphere  were  observed  over  Alaska  in 
February  of  1965  (Wong,  1965).  Thus,  it  may  be  possible  that  the  layer 
at  17  km  is  connected  with  radioactive  debris. 

The  dust  distribution  above  17  km  is  characterized  in  this  sounding  by 
layers  in  which  the  mixing  ratio  is  constant  with  altitude  but  discontinuous 
at  the  boundaries*  Another  peculiar  layer  appears  at  26  km.  The  size 
distribution  in  this  layer  is  extremely  different  from  the  surrounding 
atmosphere.  This  layer  is  also  marked  in  the  ozone  by  a  change  in  the  profile 
structure,  not  by  a  change  in  concentration. 

A  meaningful  analysis  of  the  dust  and  ozone  profile  can  only  be  made 
when  enough  data  is  available  to  identify  the  same  phenomenon  in  more  than 
one  flight.  A  comprehensive  review  and  analysis  of  the  entire  data  obtained 
from  this  projecy  is  forthcoming. 
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Solar  Extinction 
by 

T0  C.  Gillett  and  T.  J.  Pepin 

A  number  of  high  altitude  balloon  flights  have  been  made  with  an 
instrument  to  measure  the  extinction  of  sunlight  during  sunrise  as  the  path 
length  of  the  ray  in  the  atmosphere  varies  with  the  sun’s  position  on  the 
horizon  to  8°  above  the  horizontal.  These  measurements  are  similar  in 
kind  and  purpose  to  the  measurements  of  stellar  extinction  reported 
previously  (Ann.  Prog.  Report  1955).  Although  the  phenomenon  of  ’’scintillation" 
cannot  be  observed  a  measurement  of  solar  extinction  has  the  basic  advantage 
of  increased  light  level  which  allows  a  wide  choice  of  detectors  and  spectral 
ranges  for  investigation. 

These  first  experiments  were  designed  to  supplement  the  direct  measurements 
of  dust  and  ozone  and  the  extinction  was  measured  in  three  different  wave 
length  bands  which  would  best  distinguish  between  the  effects  of  these 
constituents.  Photocells  were  used  as  detectors  to  measure  the  extinction  in 
a  band  near  4500  A°,  6000  A0  and  8SOO  A0.  The  extinction  in  the  4500  A°  region 
is  primarily  Rayleigh  scattering  by  air  molecules  with  only  a  small  contribution 
due  to  scattering  by  dust.  The  band  at  6000  A0  is  at  the  maximum  of  the  Chappius 
absorption  by  ozone  and  the  extinction  in  this  band  will  be  primarily  produced 
by  ozone  absorption.  In  the  infrared  region,  the  8500  A°  band,  the  extinction 
should  be  produced  mostly  by  dust  scattering  with  only  a  small  amcunt  of  Rayleigh 
scattering.  The  data  on  solar  elevation  and  intensity  have  been  reduced  to  a 
plot  of  intensity  as  a  function  of  air  mass  and  in  this  form  (on  a  semi- log 
plot) ,  there  is  found  for  all  three  wave  length  regions  a  linear  variation  at 
small  path  lengths.  This  makes  it  possible  to  define  an  extinction  coefficient 
fcrthe  upper  atmosphere  fr^m  the  relation:  I=I0  exp  C-kX)  where  X  is  the  path 
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langth  in  units  of  air  mass.  The  extinction  coefficient  in  the  infrared 
(8500  A0  band)  was  found  to  be  highly  variable.  The  maximum  value  of 
observed  was  0,14  on  28  August  1965  and  a  minimum  of  0,07  on  11  November 
1965,  a  change  attributed  to  an  equivalent  change  in  the  number  density  of 
small  dust  particles  in  the  upper  atmosphere. 

Preliminary  analysis  has  also  been  performed  on  the  data  from  one  of 
these  flights  cn  rays  passing  through  the  atmosphere  lying  below  the  balloon 
to  determine  the  vertical  distribution  of  dust  down  to  the  tropopause.  The 
results  of  this  analysis  obtained  by  a  numerical  t rans fornat ion  of  the  data 
on  extinction  as  a  function  of  atmospheric  path  is  shown  in  Figure  1  where 
the  product  of  the  number  density  of  dust  particles  times  effective  scattering 
cross  section  for  the  dust  is  plotted  against  altitude.  That  these  values 
of  dust  concentration  are  consistent  with  the  direct  measurements  by  Rosen 
(1966)  may  be  appreciated  by  assuming  a  characteristic  size  of  Q,5w  for  the 
scattering  particles.  The  maximum  concentration  just  above  the  tropopause 
then  corresponds  to  a  number  density  of  several  dust  particles  cm-3. 

On  the  basis  of  one  calculation,  this  method  of  analysis  lodes  very 
premising  and  will  be  applied  to  the  data  collected  in  the  other  flights  of 
this  series  to  determine  the  product  of  the  dust  density  and  effective 
scattering  cross  section  as  a  function  of  altitude.  It  can  further  be  extended 
to  the  extinction  results  observed  in  the  ozone  absorption  band  to  determine 
the  ozone  concentration  as  a  function  of  altitude  and  by  proper  choice  of 
spectral  region  in  the  infrared  region  of  the  spectrum,  perhaps  the  H2O  content 
of  the  atmosphere  at  high  altitudes  can  be  measured. 
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Dust  Collection  from  the  Stratosphere 

by 

Janes  M.  Rosen 

An  effort  is  being  made  to  collect  a  small  amount  of  stratospheric 
dust  for  study  in  the  laboratory.  The  instrumentation  has  already  been 
developed  and  filter  samples  will  be  taken  in  the  near  future. 

A  comparison  of  the  results  of  this  experiment  with  those  of  the  photo¬ 
electric  dust  counter  will  shew  what  portion  of  the  particles  are  opaque 
and  perhaps  will  indicate  how  the  aerosol  concentration  in  the  stratosphere 
is  altered  by  water  vapor.  The  origin  of  the  dust  can  be  determined  by 
looking  for  cosmogenic  nuclides  induced  by  the  solar  wine.  Chemical  analysis 
however,  will  not  be  possible  with  the  small  sample  obtained. 

Figure  1  shows  the  probable  size  distribution  of  dust  at  50  mb.  This 
graph  can  be  used  to  estimate  the  largest  particle  that  will  be  collected 
on  the  filter.  Since  the  sampling  volume  is  3  x  105  cm3»  the  largest  particle 
expected  is  about  6  microns  diameter.  Note  that  a  volume  of  100  cubic  meters 
would  be  necessary  to  obtain  a  few  particles  of  20  micron  diameter. 

The  dust  on  the  filter  will  be  altered  in  two  ways  from  its  natural 
state  in  the  stratosphere:  first,  the  volatile  component  (such  as  water 
vapor)  will  have  disappeared  and  seoond,  only  the  opaque  particles  will  be 
observable , 

It  is  necessary  to  have  a  reasonable  number  of  dust  particles  per  field 
of  view  in  the  microscope  to  insure  a  reliable  count.  The  filter  area 
(1  cm2)  has  been  chosen  so  that  there  will  be  about  20  particles  per  field 
of  view  with  diameters  greater  than  ,5  microns,  Experience  has  shewn  that  the 
size  distribution  can  be  determined  reliably  down  to  .3  nr  crons  diameter  using 
standard  light  microscope  techniques,  It  is  not  necessary  to  have  as  high  a 
surface  concentration  of  large  particles  because  they  are  easier  to  observe. 
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A  schematic  diagram  of  the  sampling  instrument  is  shown  in  Figures  2 
and  3.  The  vacuun  tank  draws  air  through  the  filter  after  the  balloon 
reaches  ceiling.  The  filters  are  kept  covered  except  during  the  sampling 
period.  When  the  filters  are  covered,  the  filter  chamber  is  vented  to  the 
outside  through  a  breathing  filter.  Thus,  the  dead  air  space  behind  the 
sampling  filter  cannot  draw  in  polluted  air  during  descent.  After  sampling 
has  been  completed,  valve  #1  is  closed  and  the  vacuum  tank  is  opened  to  the 
atmosphere  through  valve  #2,  This  prevents  any  contamination  that  might 
arise  from  a  leak  in  valve  #1,  The  operation  of  the  valves  and  filter  cover 
is  monitored  and  transmitted  back  to  a  ground  receiving  station.  In  addition, 
air  temperature  and  pressure  is  measured  both  inside  and  outside  the  tank. 
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Stellar  Extinction  and  Scintillation 


Although  no  further  stellar  extinction  measurements  have  been 
made  this  past  year,  there  has  been  some  development  of  the 
observation  technique  and  analysis.  Stellar  extinction  measurements 
in  principle  may  be  used  for  the  same  purpose  as  solar  extinction 
measurements  and  have  the  advantage  of  extending  the  observing 
period  beyond  the  few  minutes  of  sunrise  or  sunset.  These 
measurements  might  be  used  to  detect;  for  example,  the  diurnal 
variation  of  a  constituent.  3ecause  of  the  much  smaller  dimens ion 
of  the  stellar  bean,  it  is  also  possible  to  observe  smaller  scale 
fluctuations  in  intensity  that  nay  be  associated  with  smaller  scale 
inhomopeneities  in  the  atmosphere.  It  is  anticipated ,  therefore, 
that  some  further  exploratory  measurements  of  stellar  extinction 


will  be  made. 
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